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Abstract 

A value of (—6.1 lb 0.9stat ± 0-^syst) X 10~^ is obtained for the quadratic slope 
parameter h in the Sir^ decay Dalitz plot at the NA48 experiment at the 

CERN SPS. The result is based on 14.7 x 10^ fully reconstructed Kl 37r° ^ 
67 decays. This is the most precise measurement of any of the Dalitz plot slope 
parameters in the charged and neutral kaon decays so far. 
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Introduction 



The X — > Stt decays Dalitz plot distributions can be expanded in powers of 
Dalitz plot variables u and v [1]: 

\M{u,v)f (xl + gu + jv + hu"^ + kv^ (1) 

(S3 - So) (Sl - S2) 
U = ^ 5 V = ^ 5 

So = ^ , Si = [Pk - Pi) , ? = 1, 2, 3 



where Pk and Pj are the four-momenta of the decaying kaon and i— th pion 
(i = 3 for the "odd charge" pion), respectively. 

In the case of the Kl — > 37r° decay the expression (1) reduces to [2]: 

(xl + hB? , (2) 
^ (^u'^ + ^;2^3^ ^ Q ^ arctan (T;/V3ii) 




because of the identical final state particles. A positive/negative value of the 
quadratic slope parameter h would mean that asymmetric/symmetric final 
states are favoured. 



Based only on isospin symmetry and some general assumptions it is possible 
to define a relation 

p = h + ?,k-—^ (3) 

4cos2/? 



between the quadratic /i, k and linear g Dalitz plot slope parameters and the 
final state interaction phases (5 which should be approximately the same for 
all K ^ ?)-n decays [3]. For — > 377*^ decays the parameter p is equal to 
2 X h {k = h/3) and is not influenced by radiative corrections since all the 
flnal state particles arc neutral. Thus a high statistics measurement of the 
quadratic slope parameter in the Kj^ — > 37r° decay Dalitz plot could establish 
strong constraints to the Dalitz plot slope parameters and f.s.i. phases in 
other K — > 37r decays. Also, it could be useful in the evaluation of radiative 
corrections applied to tt^tt^tt^ and Kl 7r"'"7r~7r° decays. On the 

other hand, by combining the Kl — > 37r° decay Dalitz plot slope parameter h 
with the linear and quadratic slope parameters in the K"^ —>■ tt^tt^tt'^ decay 
it is possible to probe the vahdity of the A7 — 1/2 rule. 
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The Kl — > 37r° Dalitz plot slope parameter h has been estimated in the 
framework of ChPT [8], [9], [10]. The analysis includes a nonzero A/ = 3/2 
amplitude in the quadratic term of the Dalitz plot expansion (1) by taking into 
account the next-to- leading 0{p'^) corrections. The phenomenological couphng 
constants in the chiral Lagrangian are evaluated using a fit to the recent 
experimental data in both X — > 27r and X — > Stt decays. 

In this paper we present a measurement based upon a data set collected with 
the NA48 detector in 1998. The detector acceptance has been estimated using 
a detailed Monte Carlo (MC) simulation. 



Detector setup and data taking 



The NA48 experiment is dedicated to the measurement of direct CP violation 
in g ^ 2ti decays. The Kl beam is produced by 450 GeV/c protons from 
the CERN SPS hitting a 2 mm diameter and 400 mm long beryllium target 
at a production angle of 2.4 mrad. Passing through a set of coUimators, a 
±0.15 mrad divergent beam enters the decay volume 126 m downstream 
of the target. The evacuated 89 m long decay volume is followed by a helium 
tank which contains four drift chambers of the charged particle spectrome- 
ter. After the tank a scintillator hodoscope, a liquid krypton electromagnetic 
calorimeter, a hadron calorimeter and muon veto counters are placed. 

The liquid krypton calorimeter (LKr) is designed to measure the energy, po- 
sition and timing of electromagnetic showers [4]. The 127 cm long detector 
consists of 13212 readout cells in a projective tower geometry which points to 
the middle of the decay volume. Each cell of 2 x 2 cm^ cross section is made 
of coppcr-bcryllium ribbons which are extended longitudinally in a ±48 mrad 
accordion structure. The readout cells are contained in a cryostat filled with 
about 20 t liquid krypton at 121 K. The initial current induced on the copper- 
beryllium electrodes is measured by 80 ns FWHM and digitized by 40 MHz 
FADCs. For the 1998 data taking period, the energy and position resolutions 
of the calorimeter were determined to be: 

aiE) (3.2 ±0.2)% (0.09 ± 0.01) , , 

-^ = - ^^e^ ^©(0.42 ±0.05)% (4) 

ti yjE E 

0.4 cm 

a\x) ~ a\y) ~ — -1=^ © 0.05 cm (5) 
V E 

with E measured in GeV. 

The time resolution was better than 300 ps for photons with energies above 
20 GeV. 



5 



The energy non-linearity in the calorimeter response was found to be less than 
1 X 10-3 between 6 GeV and 100 GeV [7]. 



The LKr detector contains the so-called neutral hodoscope which is made of a 
4 mm thick plane of scintillating fibres placed near the maximum of the e.m. 
shower development. The neutral hodoscope is used in the estimation of the 
main trigger inefficiency and for an additional measurement of the event time. 

A description of the whole experimental setup can be found elsewhere [5] . 

The Kl — > 37r° — > 67 data sample was acquired by three different triggers. 
The first one is the neutral 27r° trigger (NUT 2tt^) in which each 25 ns the 
calibrated analogue sums of 2 x 8 LKr cells are used to construct 64-channel 
X and y projections of the energy deposited in the calorimeter [6]. Based 
on these projections the total deposited energy Esum-i the first and second 
moments of the x and y energy distributions, and the time of the energy peaks 
are computed. The centre-of-gravity COG of the event and the longitudinal 
vertex position are reconstructed as: 
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CQQtrig ^ V ^g^ 

Esum 



J Esum + M2,y) - (mi, + Mly) 

Jrig _ _ _y ^ Ct^, 

-^vertex — -^LKr \ ' J 



rriK 



where Mi^, , ^h,y and M2^x^ ^2,y are the first and second moments of the 
energy deposition, ZLKr is the longitudinal position of the calorimeter and 
niK is the nominal kaon mass. 

The proper decay time r**"*^ from the beginning of the decay region (just 
after the final collimator) in units of the Ks lifetime ts is derived taking into 
account calibration constants. The events are accepted by the NUT 2n^ trigger 
if Esum > 50 GeV, COG^"<' < 15 cm and t*"V'^s < 5. In addition, less than 
6 energy peaks within 13.5 ns in both the x and y projections are required. 

The second trigger (NUT Svr^) uses the same hardware chain as NUT 2tt^ 
trigger but is specially set-up to select 37r° events by applying no condition on 
the number of the energy peaks. 

The third trigger (Nhodo) is based on the information from the neutral ho- 
doscope. It requires a coincidence of signals from the upper and lower or the 
left and right parts of the calorimeter. In order to reduce the trigger output 
rates the Nhodo trigger and the NUT 37r° one were properly down-scaled. 
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Data analysis 



The Kl 37r° — > 67 events were selected from all events which met at least 
one of the trigger requirements. The selection criteria required 6 or more LKr 
clusters satisfying the following conditions: 

• The energy of the cluster had to be between 3 GeV and 100 GeV; 

• To avoid energy losses in the LKr, the distance from the cluster to the 

closest dead calorimeter cell was required to be greater than 2 cm and the 
cluster had to be more than 5 cm from the beam pipe and the outer edge 
of the calorimeter. 

All possible combinations of 6 clusters which passed these requirements were 
considered. In addition, the following further selection criteria were applied 
on each combination: 



The distance between all cluster pairs had to be greater than 10 cm in order 
to avoid overlapping of the clusters; 

All 6 clusters had to lie within 5 ns from the average cluster time; 

The sum of the clusters energies had to exceed 60 GeV which is sufficiently 

far from the NUT 27r'' and NUT Stt^ trigger threshold of 50 GeV; 

The radial position of the center-of-gravity 



COG^ 



Ei=i Ei 



where E^, Xi, Hi are the energies and positions of the six selected clusters, 
had to be less than 10 cm; 
• No additional clusters with an energy > 1.5 GeV were allowed within ±3 
ns from the event time to minimize possible accidental effects. 

The longitudinal vertex positions were reconstructed analogous to (7): 

^ELi E5>, E,E^ [{x, - x^f + (y, - y^f] 

'^vertex — ^LKr loj 

rriK 



where and Xi, yi are the energies and positions of the six selected clusters. 

For each combination the invariant masses mi, 1712 and of all 15 possible 
photon pairings were computed using z^ertex- By applying a criteria, the 
combination and pairing most compatible with the hypothesis that mi, m2 
and rris are equal to the nominal 7r° mass were picked. 
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Table 1 

The summary of the number of accepted events and the measured values of the 
Dalitz plot slope parameter h for the NUT Stt", NUT 27r° and Nhodo triggers. The 
errors are statistical only. 





NUT 37r° 


NUT 27r0 


Nhodo 


Data events 


12.43 X 10^ 


1.48 X 10^ 


0.82 X 10^ 


MC events 


4.76 X 10^ 


4.53 X 10^ 


2.49 X 10^ 


h X 10=^ 


-6.36 ±1.06 


-4.49 ± 1.88 


-7.53 ± 2.50 



To ensure the purity of the sample, events which had one or more reconstructed 
tracks in the spectrometer or more than one hit per plane in the third and 
fourth drift chambers were rejected. 

To improve the resolution on the Dalitz plot variables R'^ and tan(^) and to 
reduce differences in the energy resolutions and non-linearities between data 
and MC which may bias the measurement of the quadratic slope parameter 
h, all the events were passed through a kinematical fitting procedure with 
constraints. Assuming the neutral kaon coming from the Kl target and a 
K — > 37r° 67 decay the cluster energies and positions are adjusted to 
minimize 

, ^ i^'-^rf ' (^.-^rf ' (y.-y!"f 

i=i cr{Ei) a{xi) a {yi) 



where a{Ei), a{xi), a (yi) arc given by (4), (5) and , xf' , yp are the 
adjusted energies and positions of the photons as functions of the 15 decay 
parameters (the kaon energy, the vertex position, the three Euler angles of 
the decay plane in the kaon rest frame, the Dalitz plot variables and 6 and 
the 2 angles of the 7 directions for each of the decaying 7r°s) . While the effect 
of the kinematical fitting on the E^ and z^ertex resolutions was small, the R^ 
and 9 resolutions were improved significantly. In particular the RMS widths 
of - i?2 and {9 - 9true) from 8.8 x lO'^ and 48 mrad became 3.0 x lO'^ 
and 31 mrad, respectively. 

Finally a cut X/jf < 8.5 was applied to exclude the tails in the x^^^ distribution 
which differ for data and MC. 

The numbers of data and MC events which passed the selection criteria are 
presented in Table 1. 

To estimate the quadratic slope parameter /i, the R^ distribution for the data 
events was corrected for the detector acceptance. For this purpose about 8x 10^ 
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Monte Carlo events were generated with h = 0. Then the quadratic slope 
parameter h was evaluated by a linear fit to the ratio of the data and MC 
distributions (Fig. 1). The fit was done separately for each of the triggers 
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Fig. 1. (a) The data (solid line) and MC (dashed line) distributions for the NUT 
Stt^ trigger, (b) A linear fit to the normalized at = ratio of these distributions. 

using the corresponding MC samples within appropriate decay regions. 

An important feature of the analysis was that the ratio of the data and MC 
distributions was fitted up to 1.9 (Fig. 1). This allowed us to significantly 
improve the stability of the result by excluding the part of the Dalitz plot 
most affected by the energy resolution and non-linearities. Since the shape of 
the detector acceptance on i?^ depends mainly on the kaon energy, the whole 
analysis was done in kaon energy bins. Combining the measured values of the 
Dalitz plot slope parameter h presented in Figure 2 and Table 1 we get a 
overall result of (-6.1 ± 0.9) x 10"^ 
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Fig. 2. The Kl Stt*^ decay Dalitz plot slope parameter h in kaon energy bins for 
the NUT 37r° (a), NUT 27r° (b) and Nhodo (c) triggers, respectively. The overall 
result and the statistical error are shown by the shaded areas. 

Systematics 



Contributions of the detector acceptance estimation, the x^^^ cut, the trigger 

conditions and the differences in the energy resolution, non-linearities and 
energy scale between the data and MC to the systematic error were considered. 

In order to study possible uncertainties in the acceptance we divided the data 
and MC samples upon the decay length t/t^, the longitudinal vertex position 

^vertex a^d the centre- of- gravity COG. For each considered trigger sample the 
maximum deviation between the average value of h from Table 1 and the 
values obtained averaging over the dependence of a second variable {t/ts, 
^vertex, COG) was taken as systematic error (Table 2). 
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Table 2 

The Dalitz plot slope parameter h in units of 10~^ for the different subdivisions of 
the data and MC. 





NUT 37r° 


NUT 27r° 


Nhodo 


In Z^ertea;bins 

X^/n.d.f 


-6.50 
14.3/9 


-4.48 
7.4/10 


-6.72 
6.9/10 


In t/ts bins 
X^/n.d.f 


-6.43 

14.6/9 


-4.58 
8.9/9 


-7.70 

7.2/9 


In COGbins 
X^/n.d.f 


-6.38 
12.4/9 


-4.64 
11.6/9 


-7.56 
5.7/9 


Ah 


0.14 


0.16 


1.0 



An important check of the reliability of the MC was done by comparing the 
data and MC Dahtz plot variable 6 (Fig. 3). A linear fit to the ratio of the 
data and MC 6 distributions clearly points to a very good agreement between 
experimental data and the simulation. 

The second considered source of systematics was the difference between the 
data and MC cluster energy resolution. Using the measured and kinematically 
adjusted cluster energies we established a strong control over the resolutions 
seen in the data and the MC. The differences in the resolutions were adjusted 
by an additional smearing to the cluster energies in the MC. The maximum 
change in h after these adjustments was found to be ±0.20 x 10~^ and was 
taken as a systematic error. A similar approach was applied in the study of 
the systematics due to the energy non-linearities leading to a ±0.17 x 10~^ 
systematic error on the h. To take into account small non-gaussian tails in the 
cluster energy resolution seen in the data, a simplified parameterization was 
used to add the tails to the MC events. The effect of the tails was found to be 
negligible. Changing the energy scale in the MC by ±1 x 10~^ we got a shift 
in h of ±0.17 x lO'^. 

The influence of the trigger conditions was studied by subdividing the data 
and MC in t/ts bins. Although that for the NUT 27r° and Nhodo trigger 
samples no systematic effect was found, the result for NUT Stt*^ trigger sample 
was affected by a small NUT inefficient region near the r*"^/Ts cut. The data 
and MC were processed excluding this region and the change in the result of 
0.1 X 10~^ gave us a systematic error due to the trigger. 

The stability of the result upon the x}it cut was checked by loosening it to 
17. No statistically significant effect on h was established. Furthermore, the 
result was found to be stable upon different cuts on the clusters energies and 
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Fig. 3. The normalized ratios of the data and MC 9 distributions for the NUT 
37r° (a), NUT 27r° (b) and Nhodo (c) triggers. The ratios were fitted with 
f{e) = {1 + Slopexe). 



positions in LKr. 

No influence of a possible background contamination from 7 conversions and 
Dalitz decays was observed. 

All the sources of systematic error were added quadratically and yield a total 
systematic error of 0.5 x 10~^. 

In addition, a second independent analysis has been performed yielding a 
consistent result. 
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Conclusions 

Our result for the Dalitz plot slope parameter in the Ki 37r° decay is h = 
(— 6.1 ±0.9stat ±0.5syst) X 10"'^, where the errors are statistical and systematic, 
respectively. This result is in an agreement with the value of (—12 ± 4) x 
10~^ predicted in [8]. The previous measurement of ( — 3.3 i \.\giat i syst) ^ 
10"'^ performed by the E731 experiment [11], is comparable with the result 
presented in this paper. 
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